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Abstract 

In many extensions of the Standard Model, including a broad class of left-right symmetric 
and Grand Unified theories, the light neutrino mass matrix is given by the left-right symmetric 
seesaw formula = fv^ — ^Yyf~^Yy, in which the right-handed neutrino mass matrix and 
the SU{2)l triplet couplings are proportional to the same matrix /. We propose a systematic 
procedure for reconstructing the 2" solutions (in the n-family case) for the matrix / as a function 
of the Dirac neutrino couplings {Yi,)ij and of the light neutrino mass parameters, which can 
be used in both analytical and numerical studies. We apply this procedure to a particular 
class of supersymmetric 5*0(10) models with two 10- dimensional and a pair of 126 © 126 
representations in the Higgs sector, and study the properties of the corresponding 8 right- 
handed neutrino spectra. Then, using the reconstructed right-handed neutrino and triplet 
parameters, we study leptogenesis and lepton flavour violation in these models, and comment 
on flavour effects in leptogenesis in the type I limit. We find that the mixed solutions where 
both the type 1 and the type II seesaw mechanisms give a significant contribution to neutrino 
masses provide new opportunities for successful leptogenesis in 50(10) GUTs. 
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1 Introduction 



Experimental data suggest that neutrinos are massive and mix. The most popular explanation 
of the smallness of their masses relies on the (type I) seesaw mechanism Qj, which finds a natural 
realization in Grand Unified Theories (GUTs) based on the 5*0(10) gauge group. While the 
seesaw mechanism cannot be directly tested since, at least in its GUT version, it involves 
superheavy states, it has observable consequences: leptogenesis j2] and, in supersymmetric 
theories, flavour 3J and CP violation |lj in the lepton sector. Successful leptogenesis puts 
several constraints on the seesaw parameters IHl El ; in particular, the mass of the lightest 
right-handed neutrino. Mi, should be larger than 0(10^ — 10^) GeV in the case of a hierarchical 
mass spectrum. It is well-known that in 5*0(10) models where the dominant contribution to 
fermion masses comes from 10-dimensional Higgs representations, the right-handed neutrino 
mass spectrum reconstructed from the type I seesaw mass formula is strongly hierarchical 
(except for special values of the light neutrino mass parameters |8j), with Mi lying below the 
minimal value for a successful leptogenesis [H]. 

However, many extensions of the Standard Model contain two sources for neutrino masses 
and leptogenesis: the type I (associated with the exchange of right-handed neutrinos) (Ij and 
type II (associated with the exchange of heavy scalar SU{2)i triplets) [IDIIII] (see also Ref. [12]) 
seesaw mechanisms. In particular, in a broad class of left-right symmetric and Grand Unified 
theories, the light neutrino mass matrix is given by the left-right symmetric seesaw formula 
My = fvL — ^^iyf~^^u, in which the right-handed neutrino mass matrix and the heavy triplet 
couplings are proportional to the same matrix /. Very often in the literature it is assumed 
that either of these two mechanisms dominates in the neutrino mass matrix. But, as discussed 
below, this corresponds to assuming specific values of the unknown seesaw parameters. In a 
general study one should encompass the situation where both contributions are sizeable and 
can be comparable in magnitude. In order to implement in an efficient way our experimental 
knowledge about neutrino masses and mixings, we need a procedure to reconstruct the matrix 
/ as a function of the Dirac neutrino couplings. This problem has been first addressed in 
Ref. where it was found that there are exactly 2" different solutions in the n-family case, 
which are connected two by two by a transformation called "seesaw duality" . 

In this paper, we use a different, more efficient procedure to reconstruct the 2" solutions, 
which employs complex orthogonal transformations and is appropriate to both numerical and 
analytic studies. For three generations of neutrinos, the eight solutions correspond to the 
different combinations of the roots of three quadratic equations. We apply this procedure 
to a particular class of supersymmetric 5*0(10) models with two 10-dimensional and a pair 
of 126 © 126 representations in the Higgs sector, and use the results to study leptogenesis 
and lepton flavour violation in these models. The spectrum of possibilities to account for the 
observed neutrino data is much richer than in the cases of type I and type II dominance, and 
the mixed solutions where both seesaw mechanisms give a signiflcant contribution to neutrino 
masses provide new opportunities for successful leptogenesis in 50(10) GUTs. 

The content of the paper is as follows. In Section |21 we present the procedure for recon- 
structing the 2" solutions for the matrix / from the light neutrino mass parameters, assuming 
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that the Dirac matrix is known, and we discuss the properties of the solutions. In section |2l 
we apply the procedure to supersymmetric 5*0(10) models with two 10-dimensional and a 
pair of 126 © 126 representations in the Higgs sector, and display the corresponding 8 right- 
handed neutrino spectra as a function of the B — L breaking scale, for various values of the 
free parameters (which include several phases). In Section EJ we compute the CP asymmetry 
in right-handed neutrino decays for the 8 solutions, and comment on flavour effects in the type 

1 limit. In Sectional we discuss the predictions for lepton flavour violating processes. Finally, 
in Section ini we give our conclusions and comment on possible extensions of the present work. 
Analytic approximations that can be useful to understand the results of the reconstruction 
procedure are given in Appendix IHl 

2 Reconstruction of the heavy neutrino mass spectrum 

In many extensions of the Standard Model based on a gauge group embedding the left-right 
symmetric group SU (2) ^ x 5'[/(2)/j x f/(l) b-l, the neutrino mass matrix is given by the following 
seesaw formula [TUl [TT] : 

= hvL - —Yjf^'Y, , (1) 

where fi and are symmetric matrices. The second term in the right-hand side of Eq. 
is the usual type I seesaw mass term, where the heavy Majorana mass matrix Mr = JrVfi is 
generated from the vev of an SU{2)r triplet A/j with couplings {fR)ij to right-handed neutrinos, 
and Md = YyV is the Dirac mass matrix (t> = 174 GeV is the vev of the SM Higgs doublet, 
to be replaced by Vu = f sin/? in supersymmetric models). The first term, known as the type 
II seesaw mass term, is generated from the exchange of an heavy SU{2)l triplet A^, with 
couplings {fL)ij to lepton doublets. The induced vev vl is related to the heavy triplet mass 
Mai^ by Vl ~ v'^vr/ M\^, which naturally explains its smallness. 

In this paper, we consider theories in which the couplings of the SU{2)l and SU{2) triplets 
are equal and the Dirac mass matrix is symmetric, so that Eq. becomes: 

= fvL - —YJ-'Y. . (2) 
Vr 

These relations arise naturally in 50(10) GUTs in which the right-handed neutrino masses are 
generated from a 126 Higgs representation jTHj (barring non-symmetric contributions to the 
Yukawa couplings, coming e.g. from a 120 Higgs representation), as well as in a broad class of 
left-right symmetric theories jTT| . 

2.1 Reconstruction procedure 

Our starting point is the left-right symmetric seesaw formula where both / and are 
complex symmetric matrices. Our goal is to determine the matrix / for a given pattern of light 
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neutrino masses and mixings, assuming that the Dirac matrix Y"^ is known in a basis in which 
the charged lepton mass matrix is diagonal For definiteness we work in the 3-family case, but 
the procedure apphes to any number of neutrino famihes. 

If Yi, is invertible, Eq. Q can be rewritten in the form: 

Z = aX- f3X-^ , (3) 

with a = Vi/my, (3 = v'^ /{m^VR) and 

Z ^ — Y-'/'MAY-'/y , X ^ Y-y^fiY^'^Y , (4) 
nil, 

where Yu^"^ is a matrix such that Y,, = yI^'^{Y},^'^')^ (in order to deal with dimensionless quan- 
tities, we have introduced a mass scale vny characteristic of the light neutrino mass spectrum). 
Since Yy is assumed to be known in the basis of charged lepton mass eigenstates, the matrix Z 
is completely determined (up to possible high-energy phases) by the choice of the light neutrino 
mass and mixing pattern. 

Assuming further that the polynomial equation det(Z — zV) = has three distinct roots 
-2^1,2,3, we can diagonalize the complex symmetric matrix Z with a complex orthogonal matrix^: 

Z = OzDiag(^i,Z2,^3)Oi , OzO^ = l, (5) 

where the Zi are complex numbers. Then Eq. ^ can be solved for X in a straightforward 
manner, by noting that X is diagonalized by the same complex orthogonal matrix as Z: 

X = OzDiag(xi,X2,X3)0| , (6) 

with the Xi being the solutions of the quadratic equation Zi = axi — (3x~^. For a given choice 
of (xi, X2, Xa), the matrix / is given by: 

/ = Y^'^X{Y^'Y = r;/2 0zDiag(a:l,X2,X3)Oi(F;/2)^ (7) 

and the right-handed neutrino masses Mj = fiVfi are obtained by diagonalizing / with a unitary 
matrix: 

/ = UffUj , / = Diag(/i,/2,/3) , UfU] = 1 , (8) 

where the are chosen to be real and positive. The matrix Uj relates the original basis for 
right-handed neutrinos, in which Y^ is symmetric, to their mass eigenstate basis. It can be 
used to express the Dirac couplings in terms of charged lepton and right-handed neutrino mass 
eigenstates, as F = f^j^y 

Since there are two possible choices for each Xj, we have 8 different solutions for the matrix 
/ (2"' in the n-generation case), a property already found in Ref. P^]- The advantage of the 
above procedure over the one presented in Ref. is that it is more systematic and allows for 
an easier reconstruction of the right-handed neutrino masses and mixings, both analytically and 
numerically. Furthermore, the connection between the different solutions is more transparent. 

^This is also the case if det^Z — zl) = has a multiple root zi, but any non-trivial complex vector v such 
that Zv = ziv satisfies v.v ^ 0. It should be stressed that, since a complex orthogonal transformation does not 
preserve the norm of states, = v* .v, Eq. Q is not a diagonalization in the physical sense. 



3 



2.2 Properties of the 8 solutions 



In order to label the 8 different solutions for /, we denote the two solutions of the equation 
Zi = axi — (3xJ^ by x'l and a;~, with 



^ Zi ± sign(Re(zi)) zj + 4a/? 

- • (9) 



With the above definition, one has, in the limit where 4a/3 <C \zi\ ■ 



xt - - , - -- , (10) 



a Zi 



while in the limit IzjP <^ 4a/5: 



±sign(Re(2i))v/^. (11) 



We label the 8 solutions for / in the following way: (+, +, +) refers to the solution {x\^ x^, X3 ), 
(+, +, — ) to the solution (x^, xj, xj), and so on. We will sometimes refer to x~ as the "type I 
branch" and to x^ as the "type II branch" , for reasons that will become clear below. 

For some particular solutions, Eq. (fTUI) has an immediate interpretation in terms of domi- 
nance of the type I or type II contribution to light neutrino masses. More precisely, for values 
of a/3 such that 4a/5 -C \z-^ , the solutions (+,+,+) and (—,—,—) practically coincide with 
the "pure" type II and type I cases, respectively. Indeed, plugging the approximate relations 
~ Zja and X^" " ") ~ -^Z'^ into Eqs. dH), one obtains: 



^(+,+,+) 4«,_<|..P (type II limit) , 



,2 



(12) 



_^_Y^M-^Y^ (type I limit) . (13) 

'Vr 

The remaining 6 solutions correspond to mixed cases where, even in the 4a/? -C \z\\^ limit, 
the light neutrino mass matrix receives significant contributions from both types of seesaw 
mechanisms. Depending on the Dirac couplings as well as on the light neutrino mass and 
mixing pattern, one may have solutions where e.g. the type I contribution dominates in some 
entries of the light neutrino mass matrix, while both seesaw mechanisms contribute to the other 
entries^. 

There is another range of values for a(5 in which / reaches a remarkable limit, namely 
12:3 P <^ 4a/5. In this region, one has xf ~ ±sign(Re(zi)) /5/a for all which indicates a strong 



■^In fact, analogously to the type I case in which the complex orthogonal matrix R introduced by Casas and 
Ibarra [Tl| in order to parametrize the seesaw mechanism can be interpreted as a "dominance matrix" |15| , one 
may define a dominance matrix Q = ^m^ [Z^^'^)'^ {yI^'^Y [M^^^"^)^ such that = Ylik Q\i^i^ where QliTUi 
is the contribution of Zk to m^. If e.g. the contribution of zi dominates in (i.e. IQisp ^ IQaaP)! 
then one can say that either the type I or the type II contribution dominates in 7713, depending on whether 
xi = Xi or xi — Xi- 
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cancellation between the type I and type II contributions to the light neutrino mass matrix. 
This can easily be seen for the two solutions labelled by (±ei, ±62, ±63), = sign(Re(zj)), in 
which one has X ~ -^^J (3/a 1. Using Eqs. (jH), this leads to: 

^(±e„±e„±.3) ±^W^Y, , (14) 

which shows that the type I and type II contributions approximately cancel in M^. For the 
other 6 solutions, Eq. (fT^ does not hold but one still has ~ \J l3/ayi (where the Ui are the 
eigenvalues of Y^^) in the \z^\'^ -C 4a/3 regime, provided that the Dirac matrix has a hierarchical 
structure. Moreover, for the type of hierarchy considered in Section El one can show that / ~ 
i/^7«f^jDiag(si?/i,S2i/2,S3l/3)f^i/ in the \z^\^ -C Aa(3 regime, where Y^ = f/jDiag (yi, 7/2, Z/3)f^i^ 
and Si = ± sign(Re(2;j)) for Appendix |B)) . 

Finally, in the intermediate region of values for ajS, {zil"^ < AajS < I-Zap, both the type 
I and the type II seesaw mechanism give significant contributions to the light neutrino mass 
matrix. Already for <^ 4a/5, cancellations between the right-handed neutrino and triplet 
contributions to Mi, start to occur. 

The 8 different solutions for / are connected to each other by the three transformations: 

which act This generalizes the "seesaw duality" of Ref. PH], defined as / 

/ = My/vL — f, which amounts to interchange the type I and type II branches for all three 
Xi simultaneously, thus dividing the 8 solutions into 4 "dual pairs". The transformations ()15p. 
more generally, allow to generate all 8 solutions from a single one. In group-theoretical terms, 
these transformations define an abelian group Z2® Zi® ^2 of 8 elements, one of which is the 
"seesaw duality" of Ref. [121 . 

3 A case study: right-handed neutrinos in 5*0(10) models 

The procedure described in Subsection 12.11 can be used to determine the a priori unknown /jj 
couplings in theories which predict the Dirac matrix Y"^, taking low-energy neutrino data as an 
input. In the following, we apply it to reconstruct the right-handed neutrino mass spectrum 
in a class of supersymmetric 5*0(10) models with two 10-dimensional and a pair of 126 © 126 
representations in the Higgs sector. 

3.1 Input parameters 



In supersymmetric SOiVS) models with two 10-dimensional and a pair of 126 © 126 represen- 
tations (but no 120-dimensional representation) in the Higgs sector, the most general Yukawa 
couplings read: 

Y^f 16il6jl0i + 16il6jl02 + 16il6jT26 , (16) 
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where Y^^\ Y^'^^ and / are complex symmetric matrices. Assuming that the SU{2)l doublet 
components of the 126 do not acquire a vev, Eq. (fT^ leads to the following mass relations for 
the charged fermions, valid at the GUT scale: 

= Md, Md = M,. (17) 

It is well-known that the second relation is in conflict with experimental data and needs to be 
corrected. In general, the corrections (coming e.g. from the SU{2)l doublet components in 
the 126 ^n]; or from non-renormalizable operators ^7j) also affect the first relation. Although 
these corrections will change numerically the solutions for /, we do not expect them to alter the 
qualitative features of our results. As a case study, we assume Eq. (fTTj) to hold in the following. 

The inputs in the procedure for reconstructing the 8 solutions for / are the matrices Y^, and 
Ml, at the seesaw scale^. The "boundary condition" for Y^, Eq. (fT7|) . is defined at the GUT 
scale, where it is convenient to work in the basis for the 16 matter representations in which Me 
(hence Md) is diagonal with real positive entries. In this basis, the Dirac matrix reads: 

Y, = U^YuU,, U, = PuVcKMPd, Yu = 'Diagiyu,yc,yt) , (18) 

where Vckm is the CKM matrix and yu,c,t are the up quark Yukawa couplings, all renormalized 
at the GUT scale. The presence of two diagonal matrices of phases Pu and Pd in Eq. ()18|1 is 
due to the fact that the 5*0(10) symmetry prevents independent rephasing of right-handed and 
left-handed quark fields. Since Y^, is only weakly renormalized between Mqut and the seesaw 
scale, the effet of the running being smaller than the uncertainty on the quark parameters at 
Mgut, we can neglect it and assume Eq. (fTHj) to hold at the seesaw scale, in the basis of charged 
lepton mass eigenstates. In the same basis, the light neutrino mass matrix generated from the 
seesaw mechanism reads: 

M, = UtM.U} , Ui = PMpmnsPu, = Diag (mi, ma, mg) , (19) 

where Upmns is the PMNS matrix and mi 2,3 are the light neutrino masses, all renormalized 
at the seesaw scale. The two relative phases in P,^ are the physical CP-violating phases associ- 
ated with the Majorana nature of the light neutrinos, while the three phases contained in Pg, 
analogous to the five independent phases contained in P„ and Pd, are pure high-energy phases. 
Having specified Eqs. (fTHj) and (fT^ . one can apply the procedure presented in Subsection 12.11 
and reconstruct the 8 different matrices / corresponding to a given light neutrino mass and 
mixing pattern as a function of a, P and of the high-energy phases contained in Pu, Pd and Pg- 

The associated right-handed neutrino mass and mixing patterns strongly depend on the 
values of a and /? (or equivalently P/a and vr), which in turn depend on the details of the 
model. The simplest way to realize the type II seesaw mechanism in the class of 5*0(10) models 
considered is to introduce a 54 representation in the Higgs sector in addition to the 126 © 126 
pair. This is easily seen in a left-right symmetric language: the 126 contains a right-handed 
triplet A'' with quantum numbers (1, 3, 1)„2 under 5f/(2)ix5f/(2)Rx5'f/(3)cxf/(l)B-L, whose 

^Actually Eq.© involves the decoupling of four states at scales that can differ by several orders of magnitude; 
the associated radiative corrections are neglected here. 
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vev vr is responsible for the breaking of B — L, as well as a left-handed triplet A = (3, 1, l)+2; 
the 54 contains a bitriplet A = (3, 3, l)o; and each 10 contains a bidoublet $ = (2, 2, l)o. The 
superpotential terms relevant for the type II contribution to neutrino masses include: 

W = ^^•LiLjA + /t$$A + AAA"A+ ... , (20) 

where the first term comes from the 16il6jl26 couplings, and the second and third terms come 
from the 10 10 54 and 54 126 126 couplings, respectively. The presence of these terms induces a 
vevvL ~ kXv'^vr/M'^^ for the SU{2)l triplet, yielding aratio (3/a = vI/{vlVr) ~ M\^/ {kXvI^. 
Depending on the other superpotential couplings, the triplet mass Ma^ may be larger or smaller 
than vr (for vr <^ Mqut, a tuning of the superpotential parameters might be necessary to 
achieve^ M^^ < vr), hence (3 /a can be larger or smaller than 1. As for vr, its value is related 
to the breaking scheme of the GUT symmetry and is also model dependent. In principle the 
B — L symmetry could be broken anywhere between the Planck scale and the weak scale; 
however the requirement of gauge coupling unification generically disfavours the breaking of 
B — L at lower scales. Detailed studies (see e.g. Ref. [1^]) have shown that the B — L symmetry 
can be broken a few orders of magnitude below the GUT scale consistently with unification. In 
the following, we allow vr to vary in the range (10^^ — 10^^) GeV. 

3.2 Right-handed neutrino spectra 

In this subsection, we display the right-handed neutrino spectra obtained from the reconstruc- 
tion of the couplings fij in the class of 50(10) models specified above. For definiteness, we 
consider the case of a hierarchical light neutrino mass spectrum with rrii = 10^^ eV, and we 
take the best fit values of Ref. [201 for the oscillation parameters (the parametrization of the 
PMNS matrix is the one adopted in the Review of Particle Properties pT]): 

Amia = m| - ml = 2.4 x lO^^ eV^ Am^^ = - mj = 7.92 x lO^^ eV^ (21) 
sin^ = 0.44 , sin^ Ou = 0.314 , sin^ = 0.009 . (22) 

The PMNS phase 6 and the two relative Majorana phases contained in Pi, are treated as 
free parameters in our study, like the high-energy phases contained in and Pe. The 

renormalization group running between low energy and the seesaw scale has little impact on 
the neutrino parameters in the case of a hierarchical spectrum, apart from an overall scaling 
of the light neutrino masses [22] which we roughly take into account by multiplying their low- 
energy values by a factor 1.2 [23- 

In the quark sector, we take into account the renormalization group running by setting 
A{Mcut) = 0.7 in the Wolfenstein parametrization of the CKM matrix, and Yu{Mgut) = 
Diag (6 X 10"^ 2.5 x 10-^ 1) x ytiMcuT), with ytiMcuT) = 0.7. In our estimate of the GUT- 
scale values for the up quark Yukawa couplings, we have taken the central values for the first and 

'"'Due to the 54126126 interaction term, the SU{2)l triplet states in the 54 and 126 representations mix 
below Vr. By suitably tuning the values of the superpotential couplings, one can decrease the mass of the 
lightest eigenvalue of the triplet mass matrix |18| . 
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second generation quark masses given in the Review of Particle Properties |2I] • Furthermore, 
we take A = 0.22, p = 0.2 and rj = 0.35, in agreement with fits of the unitary triangle 

Before presenting the 8 right-handed neutrino spectra corresponding to these inputs, let us 
mention the restrictions that apply to the reconstructed couplings fij. A first restriction comes 
from the requirement of perturbativity, i.e. the /^'s should remain in the perturbative regime 
up to the scale at which the unified gauge coupling giQ blows up, Aio ~ 2 x 10^^ GeV. As 
discussed in Appendix we can safely take /s < 1 as a perturbativity constraint at the seesaw 
scale where the fi/s are determined. One can impose a second restriction on the couplings 
fij by requiring that there be no unnatural cancellations between the type I and the type II 
contributions to neutrino masses. In practice, we shall define the fine-tuned region by: 

hsVL > F{M,),3 , (23) 

where F measures the level of fine-tuning in the (3, 3) entry of the light neutrino mass ma- 
trix: F = 10 corrresponds to a 10% fine-tuning, and so on. Such cancellations might be the 
consequence of a symmetry ensuring a proportionality relation between / and Y^, and are not 
necessarily unnatural. Nevertheless a high degree of fine-tuning would be unstable against 
radiative corrections, since this symmetry must be broken. 

We are now ready to display the 8 solutions for the couplings fij (or, more precisely, the right- 
handed neutrino masses Mi = fiVR and the entries of the unitary matrix U f that diagonalizes /) 
as a function of vr, for a given value of and of the phases 5, and $^ {i = 1, 2, 3). 

We take as a reference point the case in which the light neutrino mass spectrum is hierarchical 
with mi = 10~^ eV, [3 = a (i.e. vl = vI/vr), and the Yukawa matrices do not contain any 
CP-violating phase beyond the CKM phase. We allow vr to vary between 10^^ GeV and 10^'' 
GeV but, as we shall see, the perturbativity constraint /s < 1 generally restricts this range. 

Figs. and El show the right-handed neutrino masses Mi^2,3 and their mixing angles |(f//)i2|, 
|([//)i3| and |(f//)23| as a function of vr, for each of the 8 solutions (+, -|-, -|-) to (— , — , — ), in 
the reference case mi = 10~^ eV, P = a and 5 = = $f = $^ = $| = 0. Due to the interplay 
between the type I and type II contributions, the observed light neutrino masses and mixing 
angles are compatible with a large variety of right-handed neutrino mass spectra. As discussed 
in Subsection 123 one recovers the type I spectrum, characterized by the approximate hierarchy 
Ml : M2 : M3 ~ : : m^, in the large vr region of solution (— , — , — ). Each eigenvalue Mj 
reaches its type I value when the condition 4a/3 -C |24_ip is satisfied, i. e. when vr 3> 2f^/|z4_j|. 
This explains why Mi is constant over the considered range of values for vr, while M2 (resp. 
M3) reaches a plateau only above vr ~ 10^^ GeV (resp. f ~ 5 x 10^^ GeV). Symmetrically, the 
type II limit, characterized by the mild hierarchy Mi : M2 : M3 oc mi : m2 : m^, corresponds 
to the large vr region of solution (-I-, -|-, -f). Finally, in the small vr region {vr < 10^° GeV) 
where there is a very strong cancellation between the type I and type II contributions, one finds 
an intermediate hierarchy Mi : M2 : M3 oc mu : mc : mt (see Eq. ()14j) and discussion below). 
Although this region is not shown in the figures, one can already see that M2 : M3 ~ m^ : mt 
for Vr = 10^^ GeV. These features of the right-handed neutrino mass spectra are expected to 
hold in other models where the Dirac mass matrix has a strong hierarchical structure, but is 
not necessarily related to the up quark mass matrix like in 5*0(10) models. 
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Mi Case +++ 




10^ 
10^ 



Mi Case -++ 




10^ 
10^ 



Figure 1: Right-handed neutrino masses as a function of vr for each of the 8 solutions (+, +, +) 
to (— , — , — ) in the reference case of a hierarchical light neutrino mass spectrum with mi = 10~^ 
eV, /? = a and no CP violation beyond the CKM phase (5 = $« = $f = = $^ = 0). The 
range of variation of vr is restricted by the requirement that /a < 1. Dotted lines indicate a 
fine-tuning greater than 10% in the (3, 3) entry of the light neutrino mass matrix. 
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The differences between the 8 different spectra plotted in Fig. [T] can be understood by 
noting that each Mj can be associated (in the sense explained in Appendix ^ with one of 
the eigenvalues of the matrix X, say Xj. If Xj = xj ("type I branch"), the corresponding 
right-handed neutrino mass Mi first grows linearily with vr, then reaches a plateau for vr ^ 
2a/ a/ (3 v'^/\zj\, corresponding to xj ~ —(3/zj. If Xj = xJ ("type II branch"), Mj first grows 
linearly with vr, then grows as vj^ for vr ^ 2^ya/(3vl/\zj\, corresponding to Xj ~ Zj/a. This 
allows one to classify the 8 different solutions as follows: the 4 solutions with x^ = x^ are 
characterized by a constant value of the lightest right-handed neutrino mass. Mi ~ 10^ GeV 
(the "type I value" of Mi) over the considered range of values for vr. Among these 4 solutions, 
the 2 solutions with X2 = also have a constant value of M2 above vr ~ 10^^ GeV. The 
2 solutions with x^ = x^ and X2 = X2 are characterized by a constant value of Mi above 
vr ~ 10^3 GeV, Ml = 5 X 10^ GeV, and by a crossing of Mi and M2 close to vr = 3 x lO^^ 
GeV. Finally, the 2 solutions with x^ = x^ and X2 = are characterized by a rising Mi. 

As can be seen from Fig. ^ the perturbativity constraint forbids large values of vr in all 
solutions but (— , — , — ). The associated upper limit on vr ranges from 3 x 10^*^ to 3 x 10^^ GeV 
(except for solution (— , — , — ) for which there is no such constraint on vr), which excludes a 
breaking oi B — L above, at or just below the GUT scale. However, as discussed below, this 
conclusion strongly depends on the value of the ratio (3 /a, assumed here to be 1. If one requires 
in addition the absence of a strong cancellation between the type I and type II contributions 
to neutrino masses, the allowed values of vr are restricted to a rather small range. 

Let us now consider the patterns of right-handed neutrino mixing angles (Fig. |21). As in the 
case of the mass eigenvalues, one can recognize known limits. The type I and the type II limits 
are recovered in the large vr region of solutions (— , — , — ) and {+, -|-), respectively. The type 

I limit is characterized by small mixing angles, very close to the CKM angles, while in the type 

II limit where / — ^ M^/vl-, the right-handed neutrino mixing angles are given by the PMNS 
angles (and since ^13 is close to its present experimental upper bound in the fit that we used, 
even |(t//)i3| is relatively large in this limit). In the small vr region {vr < 10^° GeV), the mixing 
angles are close to the CKM angles in all solutions. This can be immediately understood in the 
(+, +, +) and (— , — , — ) cases, in which, as discussed in Subsection 12. 2| / tends to -^^J j3 / aYu 
when Aafi ^ l-ZsP; in the other cases, f// ~ f/J is a consequence of the hierarchical structure of 
the Dirac matrix (see Appendix |B|) . Some striking features of the mixing patterns emerge. In 
solutions (+, — , — ) and (— , — , — ), the mixing angles are almost independent of vr. In the other 
6 solutions, the mixing angles evolve from Uf ~ t/J at vr < 10^° GeV to significantly larger 
values at large vr (cancellations may occur for some specific values of vr). The (1,2) mixing 
angle is always of the order of the Cabibbo angle or greater, except in solutions (+, +, +) and 
(— , +, +) where a cancellation occurs close to vr = 2 x 10^^ GeV. 

So far we only considered the reference case of a hierarchical light neutrino mass spectrum 
with mi = 10~^ eV, (3 = a and no CP violation beyond the CKM case. It is interesting to see 
how different input parameters would affect the results of Figs. ^ and |21 In the following, we 
briefly discuss the impact on the right-handed neutrino mass spectrum of the ratio (3 /a, of the 
high energy phases and of the type of the light neutrino mass hierarchy. 
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'-'^ij Case +++ U^ij Case -++ 




Figure 2: Right-handed neutrino mixing angles as a function of vr for each of the 8 solutions 
(+, +, +) to (— , — , — ) in the reference case of a hierarchical light neutrino mass spectrum with 
mi = 10"^ eV, P = a and no CP violation beyond the CKM phase (5 = $" = $f = $^ = = 
0). The red [dark grey] curve corresponds to |(t//)i2|, the green [light grey] curve to |(?//)i3|, 
and the blue [black] curve to |(f//)23|- 
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Figure 3: Effect of (3 ^ a on the right-handed neutrino masses. The input parameters are the 
same as in Fig. [TJ except (3 /a = 0.01. 



Let us first consider the effect oi j3 ^ a. As can be seen by comparing Figs. C] and IHl taking 
P ^ a does not change the general shape of the solutions, but amounts to shift the curves 
Mi = Mi(vji) along the horizontal axis according to vji — > ^ya//3vji (an analogous statement 
can be made about the curves {Uf)ij = iU f)ij{vR)). For instance, in solution (—,—,—), the 
type I limit is reached at larger vr values for (3/a = 0.01 than for (3 = a (see the right panel of 
Fig. EI). Nevertheless the values of the Mi corresponding to a plateau do not depend on (3/a. 
In particular, in the four solutions characterized by = x^-, one has Mi ~ 10^ GeV over the 
considered range of values for vr, irrespective of the value of Finally, the value of has 
a strong impact on the allowed range of values for vr. as can be seen in the left panel of Fig. El 
the perturbativity constraint is more easily satisfied for large values of vr when /3/a <^ 1. This 
is due to the fact that the asymptotic value of I/33I in the small vr region, \J l3/a |(yj,)33|, is 
proportional to f3/a. Conversely, the case /?/a 3> 1 is excluded because the perturbativity 
constraint I/33I < 1 is never satisfied, except in the type I limit of solution (—,—,—). Therefore, 
perturbativity constrains the SU{2)l triplet mass to lie below the B — L brealing scale (which 
might require a fine-tuning in the SU{2)l triplet mass matrix for vr <^ Mgut), except in the 
type I limit. 

The effect of input CP-violating phases other than the CKM phase on the right-handed 
neutrino masses is illustrated in Fig. |3] In general, the presence of these phases only slightly 
affects the shape of the solutions, except in regions where a crossing of two mass eigenvalues 
occurs. Indeed, phases can lift isolated degeneracies between two eigenvalues (the curves repel 
one another instead of crossing), thus sensibly modifying the shape of the solution^. An example 
of this effect is shown in Fig.HJ where the solution (+, — , +) is displayed for two different choices 
of a non-zero high-energy phase, $2 — ^/^ (left panel) and $f = ir/A (right panel). These plots 
are to be compared with the corresponding plot in Fig. ^ where a crossing between Mi and 
M2 occurs at ~ 3 X 10^^ GeV. As for the right-handed neutrino mixing angles {Uf )ij, they 
are even more sensitive to input CP- violating phases than the Mj. 

^The opposite situation can also happen, i.e. input CP-violating phases can induce a crossing between two 
mass eigenvalues in cases where the corresponding curves do not intersect in the absence of phases. 
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Figure 4: Effect of high-energy phases on the right-handed neutrino masses. The input param- 
eters are the same as in Fig. ^ except $2 = (left panel), = 7r/4 (right panel). 
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Figure 5: Effect of the light neutrino mass hierarchy on the right-handed neutrino masses. The 
input parameters are the same as in Fig. Q except that the light neutrino mass hierarchy is 
inverted, with = 10~^ eV and opposite CP parities for mi and m2- 



Finally, the right-handed neutrino mass and mixing patterns also depend on the light neu- 
trino parameters that serve as an input in the reconstruction procedure, some of which are still 
unknown (mi, sign (Am^g), ^13, ^ and the two Majorana phases contained in P,^). It has already 
been shown in the type I case that particular values of these parameters can drastically modify 
the pattern of right-handed neutrino masses obtained in the generic case 0. It would be very 
interesting to investigate such effects in the case considered here; however, a general study of 
the dependence of the 8 right-handed neutrino spectra on the light neutrino mass parameters 
is beyond the scope of this paper. We just show in passing (Fig. EI) the impact of the type of 
the light neutrino mass hierarchy on right-handed neutrino masses, for the two solutions where 
the effect is the most significant. 
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4 Leptogenesis 



In the previous section, we showed on a particular 5*0(10) example that the spectrum of pos- 
sibilities to account for the experimental neutrino data in the presence of both type I and type 
II seesaw mechanisms is very rich. This has of course important implications for phenomena 
in which the presence of right-handed neutrinos and/or of a heavy SU{2)l triplet plays a role, 
such as leptogenesis and, in supersymmetric theories, lepton flavour violation. 

In this section, we show that taking into account both seesaw contributions to neutrino 
masses opens up new possibilities for successful leptogenesis in 5*0(10) GUTs. Since the model 
we consider is not fully realistic as it leads to wrong mass relations between charged fermions, 
we do not undertake a full study of leptogenesis including washout effects, but consider solely 
the value of the CP asymmetry. We do not try either to maximize the asymmetry by playing 
with all input parameters (in particular, we stick to a hierarchical light neutrino spectrum with 
mi = 10~^ eV and the best fit values (j?T|) and (j^^ for the oscillation parameters), but we 
restrict our attention to the impact of the input CP-violating phases. 

In the scenario we are considering, it is natural to assume that the lightest right-handed 
neutrino is lighter than the SU{2)l triplet. Indeed, while the perturbativity constraint discussed 
in Subsection IH. 21 requires Ma^ < vr, Mi lies several orders of magnitude below vr- Thus, one 
can safely assume that Mi <^ M/^^ , in which case the lepton asymmetry is dominantly generated 
in out-of-equilibrium decays of the lightest right-handed (s)neutrino. The CP asymmetry eAr^ = 
[r(A^i ^ IH) - r(iVi ^ IH*)] I \r{Nx IH) + r(A^i ^ m*)] receives two contributions: the 
standard type I contribution ejy^ and a contribution from a vertex diagram containing a 

virtual triplet, e^^ j2nilZZl- In the case Ml <^ M2^3 which is relevant here, they can be written 
as |27n28i: ' 



e 



nil) 



kl 



Ml 



(24) 



where Mj = —(3 Yf~^Y and M^^ = af are the type I and type II contributions to the neutrino 
mass matrix, respectively. The total CP asymmetry in A'^i decays then reads: 

, , 3 E../Im [yikyii iM,)U Ml , , 

eN. - e^^+e^^ - - ^^^^ — . (25) 

In Eqs. fj2l|l and (f25|l . the Dirac couplings are expressed in the basis of charged lepton and 
right-handed neutrino mass eigenstates, i.e. Yifc = {UjYiy)ik. Besides its obvious dependence 
on the light neutrino mass matrix and on the phases it contains, the CP asymmetry depends 
on the considered solution for the matrix / and on the input parameters (in particular on the 
phases) through their infiuence on the values of Mi and of the right-handed neutrino mixing 
angles (f//)ii {i = 1,2,3). 

The final baryon asymmetry is given by: 

— = -1.48 X lO^^r^ejVi , (26) 
s 
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£cp Case ecp Case 




Figure 6: CP asymmetry ei as a function of vr for the solution (— , — , — ) in the case of a 
hierarchical light neutrino mass spectrum with mi = 10~^ eV, /3 = a, and no CP violation 
beyond the CKM phase (left panel) or different choices of CP- violating phases (right panel). 
On the left panel, the thin lines correspond to the contribution of right-handed neutrinos (red 
[grey] curve) and of the heavy triplet (black curve). 



where 77 is an efficiency factor that takes into account the initial population of right-handed 
(s)neutrinos, the out-of-equihbrium condition for their decays, and the subsequent dilution of 
the generated lepton asymmetry by wash-out processes [Ij. For leptogenesis to be successful, 
Eq. (pUj) should reproduce the observed baryon-to-entropy ratio n^/s = (8.7±0.3) x 10~^^ [23]. 
Detailed studies of thermal leptogenesis in the type I case (see e.g. Refs. EI) have shown 
that ?7 > 0.1 over a significant portion of the parameter space; therefore thermal leptogenesis 
can succesfully generate the observed cosmological baryon asymmetry for \ ~ 10~^ (or even 
for |eiVi| ~ few x 10~^ in the case of a thermal initial population of Ni / Ni). Large efficiency 
factors can also be obtained in the presence of both type 1 and type II seesaw mechanisms |27j . 

Figs, ini to IHl show the absolute value of the CP asymmetry in A^i decays as a function of 
fij, for three representative solutions (—,—,—), (-I-, +, -|-) and (-!-,—,+). Before commenting 
on these results, let us note that an upper bound on e^^ can be derived from Eq. [771 OH]: 

1'^^' - =8^ - ^ ^ ^° [WGeVj [omeVJ ' ^^^^ 

where rrimax = max (mi, m2, 77^3). From this one can already conclude that, for a generic'' 
hierarchical light neutrino mass spectrum, the four solutions characterized by X3 = x^, which 
give Ml ~ 10^ GeV for all values of v^, fail to generate the observed baryon asymmetry from 
A^i decays (we comment at the end of this section on possible flavour effects). This is confirmed 
by Fig. IHl which shows that, depending on the values of the CP-violating phases, |eArJ ranges 
from 10~^^ to 2 X 10"^^ in the (— , — , — ) solution. As expected, the CP asymmetry is dominated 

^It has been shown in Ref. |S], in the context of the type I seesaw mechanism with a strongly hierarchical 
Dirac mass matrix, that for some special values of the light neutrino mass parameters, the right-handed neutrino 
mass matrix exhibits a pseudo-Dirac structure, making it possible to generate the observed baryon asymmetry 
through resonant leptogenesis. 
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6cp Case +++ ecp Case +++ 




Figure 7: Same as Fig. IHlbut for the solution (+,+,+). 



by the type I contribution for large values of vr, while the type I and the type II contributions 
become comparable and start cancelling each other below vr ~ 10^"^ GeV. The most noticeable 
fact here is that e^-^^ (like Mi) stays constant at its type I value even far away from the type I 
limit. 

The four solutions characterized by = , which for vr > 10^^ GeV give either Mi ~ 
5 X 10^ GeV (case X2 = x^) or Mi > 10^° GeV (case X2 = x^), look much more promising. 
Indeed, solutions (+,+,+) and (—,+,+) (case X2 = x^) yield large values of e^Vi, even in 
the absence of other sources of CP violation than the CKM phase (see Fig. I?)). However, the 
effective mass parameter rhi = (yy^)iif^/Mi, which controls the out-of-equilibrium condition 
and the wash-out due to inverse Ni decays, tends to be rather large (typically rhi ~ 10~^ eV). 
The corresponding suppression of the final baryon asymmetry can be compensated for by larger 
values of eiVi, but at the price of a heavier right-handed neutrino: one typically has le^il > 10~^ 
for Ml > 10^^ GeV. Such values of Mi are in conflict with the upper limit on the reheating 
temperature from gravitino overproduction, which depending on the gravitino mass and decay 
modes may lie between 10^ and 10^° GeV [301 • One may circumvent this problem by invoking 
a non-thermal mechanism for producing right-handed (s)neutrinos after inflation, e. g. decays 
of the inflaton field • 

Solutions (+,—,+) and (—,—,+) (case X2 = x^) are in principle better candidates for a 
successful thermal leptogenesis since they predict Mi ~ 5 x 10^ GeV, a value that can lead to a 
sufficient CP asymmetry while being marginally compatible with the gravitino constraint. As 
shown by Fig. |H1 the CP asymmetry generally reaches a plateau above vr ~ 10^'^ GeV, where 
depending on the phases it can be as large as 5 x 10~^ (interestingly enough, this may solely 
be due to low-energy CP- violating phases - see the right panel of Fig. Such values of e^Vi 
could be sufficient for generating the observed baryon asymmetry, provided that the wash-out 
processes are slow enough. However, the effective mass parameter rhi tends to be too large, 
typically rhi > 10^^ eV. Larger values of e^i can be obtained in the region where a strong 
cancellation between the type I and type II contributions to neutrino masses occur. In the left 
and right panels of Fig. |Sl the peak located at vr ^ 3 x 10^^ GeV is due to a near degeneracy 
between Mi and M2; there resonant leptogenesis [321 becomes possible. In the middle panel, the 
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Figure 8: Same as Fig. IHlbut for the solution (+, — , +), with no CP violation beyond the CKM 
phase (left panel), $2=^/4 (middle panel) and $2=^/4 (right panel). 



enhancement of e^i around vr = 10 GeV is not related to any mass degeneracy and is simply 
an effect of the phase $2- ^^is case too, the wash-out of the generated lepton asymmetry is 
strong ( 0.03 eV). 

Before closing this section, let us comment on possible flavour effects EH ESI ESI EZ] , spe- 
cializing for definiteness to the type I limit of solution (— , — , — ). The relevant quantities are the 
CP asymmetries in the decays of one right-handed neutrino flavour Ni into one charged lepton 
flavour /„, defined as ef = [T{Ni l^H) - V{Ni LH*)] / [T{Ni IH) + V{Ni IH*)] , as 
well as the parameters rhf = \Yia\'^v'^ /M^, which control the out-of-equilibrium conditions and 
the main wash-out processes. Because of the smallness of its mass, including flavour effects 
in the decays of the hghtest right-handed neutrino A^i jnHl EI| does not improve the situation; 
but it has been suggested that decays of the next-to-lightest right-handed neutrino N2 (whose 
mass is M2 ~ 2 x 10^° GeV here) might lead to successful leptogenesis, without |SH] or with 
flavour effects. One interesting possibility |SS1 is that A^2 decays generate a large asymmetry in 
a specific lepton flavour that is only mildly erased by A^^i decays and inverse decays. Whether 
this can happen or not depends on the values of the parameters 63, and which we 
give in Table ^ together with the other flavoured parameters for completeness. In the case 
considered (hierarchical light neutrino mass spectrum with mi = 10~^ eV, $2=^/4 and all 
other CP- violating phases but the CKM phase set to zero), we find that the lepton asymmetry 
is essentially generated in the tau flavour; unfortunately it is small {e^ = 1.4 x 10"^) and the 
wash-out by Ni decays turns out to be strong (m[ = 2.2 x 10~^). Different choices of the 
CP-violating phases might however improve the situation. 

The above discussion shows that taking into account both the type I and the type II seesaw 
contributions to neutrino masses opens up new possibilities for successful leptogenesis in 5*0(10) 
GUTs, even though, for the specific choice of input parameters made in this paper, the wash-out 
processes tend to be too strong. Different choices for the light neutrino mass parameters, or 
different combinations of the high-energy phases, could resolve this problem. Let us also recall 
that the results presented in this section were obtained using the mass relations ()17j) . which 
need to be corrected. The inclusion of corrections leading to realistic charged fermion mass 
matrices, e.g. from the 126 Higgs representation, is not expected to alter the gross qualitative 
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parameter / lepton flavour 


a = e 


a = fi 


a = T 


^1 


9 7x1 n~i3 


— fi X 1 0^^^ 


— 15x1 0-11 




5.6 X 10"^^ 


-1.5 X 10"^ 


1.4 X 10-'^ 


^3 


-1.8 X 10-^4 


5.0 X 10-13 


-4.5 X 10-11 


m° 


3.3 X 10-3 eV 


1.6 X 10-2 eV 


2.2 X 10-2 eV 


7712 


5.9 X 10"^ eV 


1.1 X 10-2 eV 


3.5 X 10-2 eV 




4.0 X 10"^ eV 


1.1 X 10-5 eV 


9.4 X 10-3 eV 



Table 1: Parameters that control flavour effects in leptogenesis in the type I case (large vr 
limit of solution (—,—,—)), in the case of a hierarchical light neutrino mass spectrum with 
nil = 10"3 eV, $2 = ^/4 and all other CP- violating phases but the CKM phase set to zero. 



features of the right-handed neutrino mass spectrum, but might modify the numerical values 
of Ml and of the right-handed neutrino mixing angles, hence the predictions for leptogenesis. 



5 Lepton flavour violation 

In supersymmetric extensions of the Standard Model, lepton flavour violating (LFV) pro- 
cesses such as the charged lepton radiative decays Ij —>■ lij arise from loop diagrams involv- 
ing sleptons and charginos/neutralinos. The relevant flavour- violating parameters are the off- 
diagonal entries of the slepton soft supersymmetry breaking mass matrices (m|)jj, {ml^ij and 
{jrffiLjij = A^jVd, expressed in the flavour basis deflned by the charged lepton mass eigenstates. 

If the supersymmetry breaking mechanism is flavour blind, flavour violation in the slepton 
sector arises from radiative corrections induced by the flavour-violating couplings of heavy 
states populating the theory between the Planck scale and the electroweak scale. Here we must 
deal with two kinds of such couplings®: the couplings of the right-handed neutrinos PI I40| 1^ . 
Yki (where Y = UjVu), and the couplings of the heavy SU{2)l triplet [111, fij- Integrating 
the one-loop renormalization group equations in the lowest approximation, one obtains the 
following expressions for the flavour-violating (off-diagonal) entries of the soft supersymmetry 
breaking slepton mass matrices: 

K).. ^ , (m2j,,. ^ , Al^ ^ -^^ol/e.a, , (28) 

where the coefficients Cij encapsulate the dependence on the seesaw parameters: 

C,.5:K^y,.„(^).3,//%ln(^). (29) 



^We do not consider the other sources of lepton flavour violation that can be present in supersymmetric GUTs, 
such as the contribution of colour triplets 0] or the contribution of the SU{2)fi triplet whose vev is responsible 
for right-handed neutrino masses, since they are model dependent. By contrast the right-handed neutrino and 
(assuming that M/^^ is known) the SU{2)l triplet contributions can be computed once the couplings fij have 
been reconstructed. 
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Figure 9: Coefficients C12 and C23 as a function of vr for the solutions (+, +, +) and (— , — , — ) 
in the case of a hierarchical light neutrino mass spectrum with mi = 10^'^ eV, f3 = a, and no 
CP violation beyond the CKM phase. The green [light grey] curve corresponds to IC23I, and 
the blue [black] curve to |Ci2|. The horizontal lines indicate the "experimental" constraints 
IC23I < 10 and IC12I < 0.1 (see text). 



Here Mu is the scale at which universality among soft supersymmetry breaking parameters (at 
least in the slepton and Higgs sector) is assumed. In the following, we take Mjj = 10^'' GeV, 
close to the Landau pole Aio where the theory becomes non perturbative. Neglecting the smaller 
contribution of the flavour-violating A-teim and working in the mass insertion approximation, 
one can schematically write the branching ratio for Ij — >■ /j7 as: 



BR(/,->/,7) IK).," 
oc 



BR(/j UviPj) m' 



ta^n^PFsusy, (30) 



L 



where rhj^ is the average slepton doublet mass, and Fsusy is a function of the supersymmetric 
mass parameters and of tan/5. The experimental upper limits BR (yU — > 67) < 1.2 x 10~^^ jl2j 
and BR (r — > fi'y) < 6.8 x 10~^ can then be translated into upper bounds on the Cu and 
C23 coefficients as a function of the superpartner masses and of tan f3 [33] . If we require that 
the mSUGRA parameters mo and M1/2 do not exceed ~ 1 TeV, then from Fig. 3 of Ref. 
we can read the approximate upper bounds^ IC12I ^0.1 and IC23I < 10 for a benchmark value 
of tan/5 = 10. For different values of tan/5, the upper bounds approximately scale as 10/ tan/5. 

In Fig.ini we compare the values of the coefficients IC12I and IC23I in solutions (+, +, +) and 
(—,—,—) with the "experimental constraints" IC12I < 0.1 and IC23I < 10, assuming (3 = a and 
Ma^ = Vr. The plot in the left panel of Fig. IHlis representative of all solutions but (+, — , — ) 
and (—,—,—). One finds that BR(r — /X7) lies below the experimental constraint for all 
allowed values of vr (unless tan/5 is large and/or some superpartners are light), but it could 
be accessible to future experiments for vr > 10^^ GeV. The decay /i — 67 is much closer to 
its present experimental upper limit for larger values of vr, and even exceeds it for vr > 10^^ 

^More precisely, for tan/3 = 10, one has IC12I < 0.1 (resp. IC23I < 20) for Mi < 300 GeV and 400 GeV < 
mg^ < 1 TeV if ^0 = 0, and for Mi < 500 GeV and mg^ < 1 TeV if too + M1/2, where Mi is the bino 

mass and to^j, is the average slepton singlet mass. 
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Figure 10: Effect of /5 7^ a on C12 for tfie solution (+,—,+). Tfie input parameters are tlie 
same as in Fig. IHl except tliat [3/a = l and Ma^ = vr for the left curve, while jd/a = 0.01 and 
Ma^ = O.lf/j for the right curve. 



GeV (with the assumed values of the supersymmetric parameters). Solutions (—,—,—) and 
(+,—,—) have a completely different behaviour; for these solutions both BR(r fi'-f) and 
BR (/i — s> 67) are well below their experimental uper limits for all values of vr. 

We have checked numerically that, except for the large vr regions of solutions (—,—,—) and 
(+, — , — ), the type II contribution always dominates in Cu and C23 (at least for j3 ~ a). This 
can easily be understood by noting that, due to the relation = Y^, the type I contribution is 
suppressed by the small CKM angles. The coefficients Cij thus essentially reflect the structure 
of the matrix /, and like the matrix /, they have a strong sensitivity to the CP-violating phases 
and to the ratio (3 /a. The impact of (3 /a is shown in Fig. ^| One can see that, for a fixed 
value of Vr, a lower ratio P/a results in reduced LFV rates. 

To summarize, in the considered class of 5*0 (10) models, LFV processes are dominated by 
the type II contribution in most of the parameter space of the 8 solutions. The predictions 
lie significantly below the experimental upper limits, except for the large vr region of most 
solutions, where depending on the values of the supersymmetric parameters fi ^ ej can even 
exceed its present upper limit. The expected improvement of the experimental sensitivity to 
LFV processes will strongly constrain this region. 



6 Conclusions 



The procedure presented in this paper for reconstructing the matrix / of the left-right symmetric 
seesaw mechanism from the light neutrino mass parameters can be applied in any theory with 
an underlying left-right symmetry which predicts (or at least constrains) the Dirac mass matrix. 
The 8 solutions for the fij couplings can then be used to study a number of issues in which the 
presence of right-handed neutrinos or heavy SU{2)l(r) triplets (or other heavy states embedded 
in the same GUT representation) plays a role, such as leptogenesis, lepton flavour violation, 
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electric dipole moments of charged leptons, or proton decay Some of these processes (e.g. 
/X — *• 67) put strong constraints on the 8 solutions, and might even exclude some of them. The 
reconstruction procedure can also be used as a tool to investigate the flavour structure of the 
right-handed neutrino mass matrix and to make progress in the quest for a flavour theory. 

In this paper, we applied the procedure to a particular class of supersymmetric 50(10) 
models with two 10-dimensional and a pair of 126 © 126 representations in the Higgs sector. 
We found a large variety of right-handed neutrino spectra compatible with the observed neutrino 
data, opening new possibilities for successful leptogenesis in SO (10) GUTs. We also studied 
LFV processes in these models and found large triplet contributions in most solutions, especially 
in the region of large vr values. As a byproduct of our study, we found interesting constraints 
on the breaking scale of the B — L symmetry, hence on the masses of the heavy states which play 
a role in issues such as gauge coupling unification. In particular, for /5 ~ a, the perturbativity 
constraint excludes values of vr above a few (10^"^ — 10^^) GeV, depending on the solution 
(except for the (— , — , — ) solution). LFV processes further restrict the range of allowed values 
for vr as a function of the supersymmetric parameters. Also, it is interesting to note that the 
region of values for vr that is disfavoured by fine-tuning arguments is the one in which gauge 
coupling unification is problematic. 

As mentioned earlier, cancellations between the type I and type II contributions to neutrino 
masses, rather than being accidental, could be due to some (broken) symmetry ensuring a 
proportionality relation between the right-handed neutrino and triplet couplings, / oc Fjy. This 
would be particularly relevant for some interesting possibilities for leptogenesis that occur in 
the fine-tuned region, or close to it. Interestingly enough, such a proportionality relation 
automatically follows from the embedding of the S'0( 10) model in an GUT with a 27h © 
27h and a 351' © 351 pairs in the Higgs sector, in such a way that the 351' representation 
contains the 126, 54 and IO2 representations of 5*0(10), while 27h contains lOi. With this 
embedding, both couplings 16il6jl02 and 16il6jl26 come from the 27i27j351' superpotential 
term. Assuming that the doublet that couples to up quarks in lOi does not acquire a vev, one 
obtains f cc Y^. Therefore the presence of a fine-tuning in the seesaw mass formula, rather 
than being unnatural, could point to an extended unification in Eq. 

A few simplifying assumptions were made in this paper. First, the gauge symmetry breaking 
aspects of the models (including the issues of doublet-triplet splitting and gauge coupling unifi- 
cation) were not taken into account; second, corrections to the "wrong" mass relation Md = Mg 
were neglected. There are good reasons to believe that these approximations are justified at the 
qualitative level. Indeed, the main features of the right-handed neutrino spectra are dictated 
by the strong hierarchical structure of the Dirac mass matrix, and would not be spoiled by 
small corrections to the basic 50(10) mass relations. Nevertheless a more detailed analysis 
using realistic mass relations is needed in order to obtain quantitative predictions, in particular 
for leptogenesis. Also, a more systematic scan over the input parameters (most notably the 
high-energy phases and the light neutrino masses and mixings) would probably provide useful 
information. Work along these lines is in progress. 
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A Perturbativity constraint on the fij couplings 

In this appendix, we discuss the constraint coming from the requirement that the couplings 
fij remain in the perturbative regime up to the scale at which the unified gauge coupling gio 
blows up, Aio = McuT exp [— g^%/g^y ) ] ; "where 610 is the 5*0 (10) beta function coefficient 
(this Landau pole is due to the presence of a 126 © 126 pair in the Higgs sector, which gives 
a large contribution to 610). 

Above Mgut, the running of the fij couplings is governed by the renormalization group 
equation: 

^ = I Tr(/t/)/,,. + 15 (//t/),, - ^ + ■ ■ ■ , (31) 

where the dots stand for the contribution of the other superpotential couplings, which we assume 
to play a subdominant role in the regime where the /j/s are large. Assuming a hierarchy between 
the eigenvalues, fi^2 ^ /a, one finds a critical value ff^iMcuT) = ^(95 - 26)/31 c/io(Mg[/t) 
above which /s diverges before the Landau pole Aio is reached. As an example, if the Higgs 
sector contains a 45 and a 54 representations in addition to the two lO's and to the 126© 126 
pair, one obtains Aio = 1-5 x 10^^ GeV and f 3^^^' (Mgut) — 2, where we have used aio^McuT) = 
1/24 and Mqut = 2 x 10^^ GeV. Since the running of the /j^'s below Mqut is much milder due 
to the decoupling of the heavy states, we can safely take /s < 1 as a perturbativity constraint 
at the scale where the /j^'s are determined. 

In the above, we have implicitly assumed that vr < Mqut- If this is not the case, 5*0(10) 
is broken into SU{b) at the scale vr, above Mgut] as a result the running of the /j/s above 
Mgut is slower than in the case vr < Mgut, and the Landau pole Aio is shifted towards a 
larger scale. In spite of these differences, /s < 1 remains a relevant perturbativity constraint. 

B Some useful analytical formulae 

In this appendix, we provide analytical approximations that can be useful to understand the 
results of the reconstruction procedure described in Subsection 12.11 Although we follow the 
assumptions of Section El with = in the basis of charged lepton mass eigenstates, the 
formulae presented below are more generally valid in the case of a hierarchical Dirac matrix Yy^ 
with eigenvalues 2/1 <^ ?/2 ^ 2/3 and small mixing angles. 

Let us first perform the diagonalization of the matrix Z, assuming for definiteness that 
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the light neutrino mass spectrum is hierarchical (mi < m2 <^ fn^). Using the notations of 



Subsection 13.11 and choosing Y^^'^ = UTYu"^., we have: 



9 



Z = Y-'I^NY-'I' , N ^ {U;UtMMlul)/ms . (32) 

Since m2/m^ ^ 0.2 and two of the three lepton mixing angles are large (with S23 ~ 

S12 ~ and Si3 < 0.2), the matrix N has a very moderate hierarchy. We can parametrize 

it as: 

(a h c \ 
h d e \ , (33) 
c e /; 

with |a|, |c| ~ m2/m3, \d\, |e|, |/| ~ 1 and \df — e^\ ~ 1712/^13- It is convenient to define the 
quantities: 

An = df-e\ A^2 = bf - ce , A13 = be - cd , (34) 
A22 = af - , A23 = ae-bc , A33 = ad - b"^ . (35) 

All Aij are of order m2/m3, and |det = |aAii — 6A12 + CA13I = mim2/m\ < {rn2/m^Y. 
The hierarchical structure of the matrix Z is essentially determined by the up quark Yukawa 
couplings: 

a/Vu b/y/y~y; c/ y/y~y't \ 
b/y/lh^c d/vc e/y/y^t . (36) 
c/^/y^t e/^/y~irt f/yt J 

The roots of the polynomial equation det(Z — zl) = are obviously all distinct, hence Z can 
be diagonalized by a complex orthogonal matrix Oz'- 



Oz \ Z2 \ Ol , (37) 




detiV _i A33 _i _i 

Zl = — r yt , Z2 = y^ , Z3 = ay^ , {3i 

A33 a 



(39) 



\ 1 ft Vvc/yt f Vvu/yt J 

where we have ordered the Zi in such a way that \zi\ < 1^2 1 < I ^3 1- In Eqs. (jHHj) and (jH^ . the 
neglected terms are of relative order yu/yc, yc/yt with respect to the dominant terms. 

We can now reconstruct the 8 solutions for the matrix /. For a given choice of {xi, X2, 
X3), the matrix / is given by / = U'^Yu^'^XYu^'^Uq, with X = O^Diag (xi, X2, X3)0|^. The 
eigenvalues fi are obtained by diagonalizing / with a unitary matrix Uf, Eq. (jSJ. Alternatively, 
one can diagonalize the matrix X = Yu^'^XYu^'^, which is related to / by / = UjXUg, and has 
therefore the same eigenvalues: 

X = U^fUl , U^U\ = 1 , (40) 
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the unitary matrix that brings / to its diagonal form being given hy Uf = UjUj^. It does 
not seem to be possible to derive simple analytical formulae for the /j that would hold for any 
value of a and (3. However, one can easily obtain approximate formulae in the regions of {a, (3) 
values where the Xi satisfy some hierarchy requirements, as we show below. Let us first define 
the following quantities: 

xi = ytXi , X2 = Vc X2 , X3 = yu X3 . (41) 

The are given by the roots of the characteristic equation: 

det{XX^-f!l) = -ff + Cft-C^f^ + C^ = 0, (42) 

with (up to subdominant terms of order yu/Uci VclVt in the coefficients of the Xi monomials): 

C\ = r-p5 F3 + — ^5 F2 + \xi\ 



+ ^i^+^j '^'^^[^) ^'^-^n^j "■"'^"'^ 

- - |2 , / I"! n \ I- - |2 , I- - |2 



lAisI 




IA23I 


\' + \ 


IA33I 


to 




IA33I 


2 




C2 = [ ^ — ) 1x2X3] + 1^ ^ — J IX1X3I + \xiX2\ (43) 

xix^lxsl +(^;r) X1X3IX2I +(^) X2a;3|xi| + c.c. 

Using the fact that = fl + /| + /|, C2 = /f /| + /f /| + /|/| and = /f /|/|, one 
immediately sees that, when there is a significant hierarchy between the Xj, the are given by 
the Xi times an order one coefficient. More precisely, one has: 

fl ^ \/CVC^ ~ ^min , /2 ^ A/C'2/Cl ~ X^iddie , /3 - \/C^ ~ ^max , (44) 

where Xmax = max(|xi|, |x2|, jxsl), Xmi„ = min(|xi|, |x2|, |x3|) and x^iddie is the remaining Xi. 
In the case -C |x2| <^ |a;3|, this reads: 



f IA33I" 1^ I f ^ |ap(|Ai3p+|A23p + |A33|^) 1^ 

- |Al3|2 + |A23p + |A33|2 l-^H ' ^2 - 2 + |b|2 + |c|2) | A33 12 1^2 

h ^ ' J2 ' 1X31 , (45) 



while in the case |x3| <^ \x2\ <^ |xi|, one has: 

h - l^il , /2 - IX2I , fl ^ Ixsl , (46) 



etc. When the hierarchy between the Xi is not so pronounced, Eq. (j44|) is no longer a good 
approximation. In the case Xmin ^ ^middle < Xmax, however, one still has: 

fl ^ \fC^2 ~ Xrmn , (47) 
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while in the case Xmin < x^iddie < x 

max • 



Ci ~ Xr, 



(4J 



The formulae for the mixing angles are more involved, but they simplify for some Xi hierar- 
chies. It is convenient to write X as: 



X 



1 


b 






/ 


Hal 


_b Vu 


h 

a 


0? 




+ X2 




_b Vu 
a yc 


1 


c 


be 


t ) 




\ - 


bA23 yu 


A23 


a 


a? 


a? / 




aA33 ya 


A33 



+ Xi 



( 



\ 



■^13 



yi 



A13A23 yuyc. 



, A13 A23 



A^ 

.33 
Al3 

A33 



yu. 
yt 



VuVc 

~yf 



4, yl 
yl 

A23 yc_ 

A33 yt 



yt 



bA 



23 



aA33 

A23 
A33 



yu_ \ 
yc * 



^23 



■^33 



yu \ 



A13 Vu 

A33 yt 
A23 'yc 
A 33 yt 
1 



(49) 



where subdominant terms of relative order yu/ Uc and yd yt are understood in each entry of the 
matrices multiplying xi, X2 and X3. Let us first consider the hierarchy \x2\ <^ |x3|. In this 
case, X is dominated by the contribution proportional to 0:3 in Eq. (I49j) . and f^x is given by: 




''12 
"12 




(50) 



where {a,b,c) = {a,b, c)/ ^y\a^~+~\b^~+~\c\^, a = ^/T" 



laP, 



P is a diagonal matrix of phases. 



|ai2p+ |^i2p = 1, and ai2 and 612 depend on the subdominant terms in X. After multiplication 



by Ug, this gives: 



KUfU 

\{Uf)23\ 



i2 + (c^ai2-a^6&i2)K,e^(*^-*")/a| > A 



lab 



(51) 
(52) 
(53) 



Thus, the hierarchy |a;i|, \x2\ -C {xsl leads to large right-handed neutrino mixing angles (can- 
cellations are possible in Eqs. though). 

Let us then consider the hierarchy {x^l <^ \x2\- In this case, depending on the value 
of X3/X2, X is dominated either by the contribution proportional to X2 alone, or by both 
contributions proportional to X2 and X3. Ux is then given by: 



— 



1 

-&13A33 

"^13^23 







A 



23 

A 



33 




^13 



bA33 yu 

1 yc _ 
(cA^3+6Ag3)A|3 X3 



n2 



>12 







-b\2 ai2 P 










1 




< 


yu 
yc 


X2 




> 


yu 
yc 


1^2! 



(54) 



(55) 
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where (A23, A33) = (A23, A33)/ -^|A23p + IA33P, and ai2 and 612 depend on the subdominant 



terms in X. After multiphcation by , this gives: 



\iUfU 

mf)23\ 



\bi2 + ai2 AjgKdC 

|A33Krf + 6i3e-*(*^ 
IA33I ~ 1 . 



I ^ ^ 
> A , 



(56) 
(57) 
(58) 



Thus, the hierarchy |x3| -C |x2| leads to large right-handed neutrino mixing angles as 
well, but cancellations are possible in (t//)i2 and (t//)i3. The same conclusion holds in the 
qualitatively similar cases <^ |a;3| < \x2\ and \x3\ -C < \x2\. 

Let us now turn to the case \x2\, 1x3 1 -C which contrary to the previous ones leads 
to small mixing angles. In the case of a strong "inverted" hierarchy 1x3 1 -C |x2| -C |xi|, with 
1^3 1 < iVu/yt) 1^1 1 and IX2I < {Vc/yt) |^i|, one obtains: 



Ux 



1 

_b*_ £3 _j_ ^1 1" 



6* A 



a* AI 



23 ^3 Vc 
yt 



a y<: 
+ 



yt 



b X3 b Vu 

a X2 a yc 

1 

^23 yc 
A*3 yt 



Ai3 Vu \ 

A33 yt » 

A23 yc 

A 33 j/t 



1 



(59) 



U^. Hence the right-handed neutrino mixing 



After multiplication by Uj, this gives Uf 
angles are given by the CKM angles, up to corrections of order yu/Vc, Vc/vt and X3/X2 (the 
same conclusion holds for |x3| < {yu/vt) |xi|, |x2| < {Vc/yt) Explicitly, one has: 

(60) 
(61) 
(62) 

SsVu, S2yc, siyt) (st = ± sign(Re(2;i)) 



I(t^/)l2| ^ 




/6 X3 


_ b yu\ 






\a X2 


a Vc) 








Vtd- 


A23 yc 
A33 yt 








HUM ^ 


Vts- 


A23 yc 
A33 yt 









The limit {z^]'^ <^ Aa(3, which yields \/a/(3 (x3,X2,Xi 
for xf), deserves a particular discussion. In this case: 



U- ~ 
— 



v 



' a* yc 



'£3 
^^2 



1 



1)4 + ^(^-1)^ 

^ S2 'a* S2 ^ SI ' AS, 



yu_ 

yt 



byu 
a yc 

^23 yc 

yt 



33-S2feA23+aiaAi3 yu 

siaAas yt 
S2 — \) A23 Vc 
Sl ' A33 yt 



(63) 



For the two solutions characterized by Si = S2 = S3, cancellations occur in the off-diagonal 
entries of Ux, implying Ux — 1 and / ~ si^J 13 / aYu, consistently with Eq. ((T3j). For 
the other six solutions, one still has Uf ~ U^ (but with larger corrections), hence / ~ 
f/jDiag (si^/i, S2I/2, s^y^)Uq. Finally, when the hierarchy is milder, with e.g. |x3| ^ 
(yu/yt) |xi| or |x2| ^ iyc/yt) |xi|, or both, Uf deviates more significantly from U^ and is 
characterized by larger mixing angles (but cancellations may occur for specific values of the 
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ratios Xi/xj). For instance, when \x^\ < {iiu/yt)\x\\ and |VisXi| ^ \x2\ ^ |a;i|, one finds 
|([//)i2| ^ iKdl, |(f//)i3| - |Kd(A23/A33)a;2/a;i| and |(f//)23| ^ |(A23/A33) a;2/a;i|. 

With the above formulae, it is possible to explain most features of Figs. Q and El For 
example, the fact that Mi takes a constant value over the considered range of values for vr in 
all four solutions with x^ = just follows from Eq. ()47|) . As can be easily checked indeed, the 
conditions 4a/5 <^ {z^]'^ and \x^\ <^ \xf\ are always satisfied for (5 javR ^ 10"'^° GeV, 
implying Mi ~ |a;j|f_R — m'^{MGUT) / i^'f^s) , namely Mi ^ 10^ GeV for a ~ 0.2. Similarly, in 
the 2 solutions with xs = X3 and X2 = X2 , the conditions 4a/5 <C |z2p and 1x2 | <^ |xf |, IX3 | are 
satisfied for ^/JJolvr > 10^^ GeV, implying M2 ~ \x2\vR ~ ml{MGUT)/mz ~ 2 x 10^ GeV. As 
for the right-handed neutrino mixing angles, the fact that they are approximately independent 
of vr and very close to the CKM angles in solutions (+, — , — ) and (— , — , — ) is due to the 
hierarchy |x3| < {yu/Vc) \x2\ < {Vu/yt) which holds over the considered range of values 
for Vr. For large values of vr, both solutions have a strong hierarchy of the Xj's, |x3| <^ 
{jju/yc) \x2\ <S {yu/yt) and the {Uf)ij are given by Eqs. (jHIH) to (jHSI)- In the other 6 
solutions, one recovers \{Uf)ij\ ~ \(ycKM)ji\ in the small vr region (corresponding to the limit 
|-Z3p <C 4a/3, namely \/ (ijavR ^ 10^'^ GeV), while the large vr region, where the hierarchy 
12^215 |2;3 1 <^ |xi| is no longer satisfied, is characterized by larger values of the right-handed 
neutrino mixing angles. 
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